Type IIP Supernova 2009kf: Explosion Driven by Black Hole 

Accretion? 
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ABSTRACT 

Unusually bright type IIP supernova (SN) 2009kf is studied employing the hy- 
drodynamic modelling. We derived optimal values of the ejecta mass of 28.1 M , 
explosion energy of 2.2 x 10 52 erg, and presupernova radius of 2 x 10 3 R & assum- 
ing that 56 Ni mass is equal to the upper limit of 0.4 M . We analyzed effects of 
the uncertainties in the extinction and 56 Ni mass and concluded that both the 
ejecta mass and explosion energy cannot be significantly reduced compared with 
the optimal values. The huge explosion energy of SN 2009kf indicates that the 
explosion is caused by the same mechanism which operates in energetic SNe Ibc 
(hypernovae), i.e., via a rapid disk accretion onto black hole. The ejecta mass 
combined with the black hole mass and the mass lost by stellar wind yields the 
progenitor mass of about 36 M . We propose a scenario in which massive binary 
evolution might result in the SN 2009kf event. 

Subject headings: stars: evolution — stars: massive — supernovae: general - 
supernovae: individual (SN 2009kf) — gamma-ray burst: general 

1. Introduction 

Type IIP supernovae (SNe IIP) are believed to originate from 9 — 25 M stars which 
end their life with the core collapse i nto neutron star a nd the subsequent envelope ejection 



via the neutrino-drive n mechanism (IHeger et al.ll2003l ). or alternatively, via the magneto 



rotational mechanism (IMoiseenko et al.ll2006l ). In fact, limits for the mass range of SNe IIP 
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progenitors (main-sequence stars) are fuzzy, because the theory of both stellar evolution 
and SN IIP explosion is still in progress and cannot predict reliably the outcome. On the 
other hand, the observational constraints for the progenitor masses are rather ambiguous. 
Indeed, the progenitors of SNe IIP reco vered from pr e-explosion images turn out to be 
predominantly low-mass (8 — 17 M ) stars (jSmarttll2009l ). while the hydrodynamic modelling 
of the handful of well-observed S Ne IIP indicates that SN e IIP primarily originate from high- 
mass (15 — 25 M Q ) progenitors (jUtrobin fc Chugaill2009l . and references there). 



More massive (25 — 100 M Q ) stars en d up with the blac k hole formation which may be 
accompanied with a yery weak explosion ( IHeger et al.ll2003l ); they are associated with un- 
derluminous SNe IIP (ITuratto et al.lll998l ). Alternati yely, the underluminous SNe IIP might 



originate from the low- mass end of massive star rang e (jChugai fc Utrobinll2000l ; iKitaura et al. 
20061 : llJtrobin et aD l2007t llJtrobin fc Chugail l2008f ) . A tiny fraction (order of 1%) of stars 



19991 : MacFadven et all 12001 



producing black holes is responsible for energetic SNe (hypernovae) powered by the disk 
accretion onto black hole (GalamaetaL 1998 



wamoto et al.lll998l : iMacFadyen fc Woosley 



Heger et al.ll2003l ). It should be stressed that all the hyper- 



novae so far observed are exploding Wolf-Rayet (WR) stars lacking any hydrogen in the 
outer layers. 

Recently a dist ant (z = 0.182) ty pe IIP SN 2009kf has been discovered by the Pan- 



STARRS 1 survey (jYoung et al.l 120091 ) which turns out to be unusually luminous for this 
class of SNe. With the mid-plateau magnitude My —18.4 mag it is 1.5 — 2 mag brighter 
compared to normal SNe IIP. Despite a larg e distance this SN be came a subject of the 
detailed photometric and spectroscopic study ( jBotticella et al.ll2010l ). Apart from the high 
optical luminosity the SN is unusually bright in the ultraviolet (UV) band of GALEX with 
the near UV absolute magnitude M NT tv = —21.5 mag a nd has unusually large photospheric 
velocity ~ 9000 kms -1 on day 61 (jBotticella et al.ll2010l ). Authors suggest that the standard 



model of SN IIP is not applicable to SN 2009kf. They mention several possibilities to 
account for this unusual SN IIP: a huge explosion energy (> 10 52 erg), a large pre-SN radius 
(> lOOOi?©), a large 56 Ni mass, and a strong interaction of ejecta with a dense circumstellar 
(CS) shell. 

Given unusual characteristics of SN 2009kf we present results of the hydrodynamic 
simulations of this object in a framework of the explosion of a red supergiant (RSG). We 
recover basic parameters of the event and find that the required explosion energy is indeed 
tremendous, comparable with the explosion energy of hypernovae. We evaluate the role of 
alternative power sources, viz. radioactive 56 Ni decay and CS interaction, and find them 
irrelevant. Implications of our results for the pre-SN and the progenitor are discussed. 
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2. Hydrodynamic Model and Presupernova 



The spherically-symmetric hydrodynamic code with one-group radi ation transfer i s em- 
ployed to model SN 2009kf. The code was described in details earlier (lUtrobinl 12004) and 

and ref- 



2009 



was used to study several other well-observed SNe IIP (lUtrobin fc Chugail 
erences there). The explosion energy is modelled by the supersonic piston applied at the 
mass-cut which presumably is a border between collapsing core and ejected mass. A non- 
evolutionary massive RSG in hydrostatic equilibrium is used as a pre-SN model. The helium 
core of pre-SN is presumably mixed with the hydrogen envelope and the density jump be- 
tween the helium core and hydrogen envelope is essentially smoothed. Arguments in favor 
of a non-evolutionary pr e-SN for the hydrodynamic models of SNe IIP were presented by 
Utrobin fc Chugail hoO$ \. 



Principal fitting parameters, which determine the light curve and the expansion velocity 
at the photosphere, are the explosion energy, ejecta mass, pre-SN radius, and amount of 56 Ni. 
" Second order" but also important ingredient of the model is a mixing between helium core 
and hydrogen envelope. The mixing determines the shape of the plateau at the fin al stage, 
and g enerally a high degree of mixing is needed to account for the light curve shape (lUtrobin 
20071 ). A mixing of 56 Ni also affects the light curve at the transition from the plateau to the 
radioactive tail. The issue of the depend ence of the hydrodynamic model on parameters was 
explored in detail earlier ( lUtrobin! 120071 ). 



Preliminary computations of SN 2009kf light curves and expansion velocities showed 
that the explosion energy in the optimal model should exceed 10 52 erg. That enormous 
energy is beyond the capabilities of the explosion mechanisms usually associated with the 
core collapse in to the neutron star and seems to indicate the same explosion mechanis m as 
for hypernovae ( IMacFadven fc Wooslevlll999l : iMacFadven et al.ll200ll ; Irleger et al.ll2003l ) . We 
therefore consider below models with the large mass of collapsing core that forms presumably 
a black hole with the rest mass of 4.5 M . The latter is close to the minimal value which 
provides the required explosion energy by the accretion of ~ 2 M & for ~ 1% efficiency. 



3. Results 

The extensive study of a parameter space led us to the optimal model of SN 2009kf 
that can describe simultaneously the V light curve (Fig. [1]), the early strong bolometric 
luminosity peak along with the end phase of the plateau (Fig. [2]), and the ex pansion velocity 
at the photosphere (Fig. [3]). The observational V light curve is taken from iBotticella et al. 



(120 lOl ). and the bolometric light curve is constructed with a black-body fit to the g, r, i, z 
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photometry and the NUV flux, only for the first three epochs, assuming a total extinction 
Ay = 1 mag. The v elocity at the photosp here of 7350 ± 350 kms -1 is recovered from the 
spectrum on day 61 ( Botticella et al.ll2010l ) using the modelling of the Ha and He I 5876 A 
line profiles. In fitting the light curves we rely primarily on the V band data rather than on 
the bolometric light curve in which the initial peak is mainly determined by the NUV flux 
which in turn is severely hampered by errors in the total extinction. 



The parameter set of the optimal model consists of the ejected mass M env = 28.1 M e 



0; 



explosion energy E = 2.2 x 10 52 erg, pre-SN radius fl n = 2 x 10 3 and the 56 Ni mass of 



0.4 Mq. The latter is the upper limit reported by iBotticella et al.l ( j2010l ). On the basis of 
the V light curve and the velocity at the photosphere we determine the SN parameters with 
accuracy of ~ 10%. The ejecta mass combined with the adopted mass of the black hole of 
4.5 M gives a pre-S N mass of 32.6 ± 3 m aximal among other pre-SNe for previously 
modelled SNe IIP (cf. lutrobin fc Chugail . 



It should be emphasized that for the adopted value of the 56 Ni mass the ejecta mass and 
the explosion energy of the model are minimal; a decrease of the 56 Ni mass would require 
a larger ejecta mass and explosion energy. The V light curve for the model with the 56 Ni 
mass of 0.0765 M & and the model without 56 Ni illustrate this point (Fig. [TJ): with the lower 
56 Ni mass the plateau is shorter, so a larger ejecta mass and explosion energy are needed 
to account for the plateau duration for the same photospheric velocity. Computations show 
that in order to produce sensible fit of the bolometric light curve with 0.0765 M Q of 56 Ni 
one needs to increase the ejected mass and explosion energy by a factor of ~ 1.5. The light 
curve at the end of the plateau depends also on the extent of the 56 Ni mixing. Remarkably, 
in the model with the 56 Ni mass of 0.4 M Q the 56 Ni should be mixed homogeneously in the 
inner ejecta up to 7700 kms -1 . 

The luminous broad initial peak of the bolometric light curve of SN 2009kf is related, 
as usually for SNe IIP, with the large pre-SN radius. The recovered radius Rq ~ 2 x 
10 3 Rq is larger than the radii of pre-SNe of well-studied SNe IIP, which lie in the range 
of 35 — 1500 R e (Table CD). Keeping in m ind that the radius of massive RSG increases 
with the progenitor mass (IHeger et al.lll997l ). the large pre-SN radius can be considered as 
an independent evidence in favor of the relatively large pre-SN mass compared to normal 
SNe IIP. 

The primary uncertainty of the derived SN parameters is related with the poorly deter- 
mined extinction in the host galaxy. We explore this uncertainty assuming a total extinction 
to be Ay = and Ay = 1.5 mag. In the case of Ay = mag the upper limit of 56 Ni mass 
is 0.16 M & . With this amount of 56 Ni the pre-SN radius should be smaller, Rq = 200 Rq, 
while the ejecta mass and the explosion energy of the model increase up to ~ 40 M & and 
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~ 3 x 10 52 erg, respectively. Note that the E/M ratio is to be invariant because the velocity 
at the photosphere on day 61 is fixed. For Ay = 1.5 mag the upper limit of the 56 Ni mass 
is 0.63 Mq. This amount of 56 Ni permits us to obtain a satisfactory fit of the light curves 
to the observations with the same ejecta mass and the explosion energy as for the optimal 
model, but with the pre-SN radius increased by a factor of ~ 2. However, if the adopted 
56 Ni mass is lower than the upper limit, then one needs to increase the ejecta mass and the 
explosion energy to reproduce the light curves together with the velocity at the photosphere 
on day 61. We thus conclude that the ejecta mass and explosion energy cannot be lower 
than the values for the optimal hydrodynamic model. 



4. Discussion 

The major result of the hydrodynamic modelling of SN 2009kf is the conclusion that the 
SN has been caused by the energetic explosion of a massive and extended RSG. Parameters 
of hydrodynamic model for SN 2009kf are listed in Ta ble [J together with parameters of other 



SNe IIP studied earlier (see lUtrobin &: Chugaill2009i ). The columns in the order contain SN 
name, pre-SN radius, ejecta mass, explosion energy, 56 Ni mass, maximal velocity of the 56 Ni 
mixing zone, minimal velocity of the hydrogen matter, mass of the collapsed core, pre-SN 
mass, mass lost by the wind, and progenitor mass. In the case of SN 2009kf progenitor 
we suggest a single star scenario, while the total mass lost through the wind outflow is 
adopted to be equal to that lost by the progenitor of SN 2004et. SN 2009kf stands out from 
the list of other SNe IIP by its extremely large explosion energy which is 10 times greater 
than the explosion energy of SN 2004et, the most energetic event among previously studied 
SNe IIP. An exceptional nature of SN 2009kf is emphasized by its position on diagrams of 
the explosion energy and the 56 Ni mass vs. progenitor mass (Fig. H]). Despite the fact that 
SN 2009kf does not deviate from the general trends on both plots the explosion physics of 
SN 2009kf is likely essentially different from that of other SNe IIP. 

Before the further discussion of implications of the results we should consider alternative 
explanations for the unusually high luminosity of SN 2009kf. Given the rather confident 
upper limit for the 56 Ni amount, the only remaining possibility for the power source is the 
CS interaction. This mechanism suggests that the kinetic luminosity, released in the radiative 
shock wave, is emitted by the cool dense shell (CDS) at the ejecta/CS interface likewise it 
is the case for SNe Iln. Quick look at the spectra of SN 2009kf at the plateau stage shows 
that the strong CS interaction is unlikely in this case, because the spectra have little to 
do with the SN Hn spect ra. Indeed, the latter never show pronounced broad absorption 



lines (cf. iFilippenkol 119971 ) ; by contrast, the Ha line of SN 2009kf has the deep absorption 
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component (see iBotticella et al.ll2010l ). quite similar to other SNe IIP. To strengthen further 
arguments against the CS interaction we note that the maximal velocity detected in the 
Ha absorption component of SN 2009kf on day 61 is ~ 13000 kms -1 . This means that 
the double-shock structure, formed by the forward and reverse shocks with the CDS at the 
contact surface, expands with the velocity > 10 4 kms -1 . Meanwhile, the Ha absorption 
component indicates that the continuum source resides at velocities v < 6000 kms -1 , which 
is in odd with the assumption that the bulk radiation originates from the high-velocity 
ejecta/CS interface. This evidence thus rules out the CS interaction as a major source of 
the SN 2009kf luminosity. 

Regardless of the role of the CS interaction, the upper limit of the wind density could 
be valuable. The wind density can be estimated usi ng the upper limit of the luminosity 
at the nebular stage Mr > —16.3 mag on day 236 ( IBotticella et al.ll2010l ). which can be 



converted into the upper limit of the optical luminosity, L < 5 x 10 41 ergs -1 , assuming 
that the spectral energy distribution is similar to that at the end of the plateau stage. To 
recover the upper limit of the wind density w e employ the interaction model based on the 
thin shell approximation ( jChugai et al.l 120071 ) for the ejecta mass and the explosion energy 
found for SN 2009kf. For the adopted 56 Ni mass of 0.4 M & the upper limit of the wind 
density parameter turns out to be w = M/u = 6x 10 15 gem -1 . This estimate refers to 
the epoch of 236 days and the radius of 6 x 10 16 cm. Assuming the wind velocity u — 10 
kms -1 we conclude that the upper limit of the wind density corresponds to the mass-loss 
rate M < 9 x 10 -5 M Q yr -1 at the stage ~2x 10 3 yr before the SN explosion. 

The question arises on the nature of the central " infernal machine" responsible for the 
huge explosion energy E fh 2 x 10 52 erg, one order magnitude exceeding the energy of normal 
SNe IIP. In fact, the explosion energy of SN 2009kf is comparable with the explosion en- 
ergy of energetic SNe Ic (hypernovae), associated with gam ma-ray bursts (e.g. SN 1998bw). 
According to the recent compilation of iTanaka et al.l (120091 ). the explosion energies of hyper- 
novae lie in the range of (0.6 — 5) x 10 52 erg; SN 2009kf falls into this interval. The absence 
of alternative possibilities compels us to suggest that SN 2009kf is caused by the "engine" 
associated with the hypernova phenomenon. The possibility that the collapsar mechanism 
could operate in mas sive stars with retained h ydro gen envelope and prod uce a very powerful 
SN was proposed bv iMacFadven et al.l ( I200lh and IWooslev et al.l feooj ). 



The widely shared vi ew is that the hypernova explosion is powered by a rap i d disk accre 
tion into the black hole ( IMacFadven fc Woosleyl Il999l ; IMacFadven et al.l l200ll ; iHeger et al. 
20031 ). Two conceivable scenarios are proposed how the black hole and accretion disk system 
with the high accr etion rate migh t arise: (i) a core collapse in a single rapidly rotating star 
(collapsar model) ( iWoosleyl Il993l ; IMacFadven fc Woosleyl Il999[ ). or (ii) a black hole (neu- 
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tron star) merger with the He-core of massive companion (merger model) (jFryer fc Woosley 



19981 ). The merger scenario suggests that the neutron star, merged with He-core as a result 



of a super-Eddington accretion, rapidly grows a black hole surrounded by the accretion disk. 

The presence of a massive hydrogen envelope in SN 2009kf imposes constraint on the 
single star scenario: the progenitor should be less massive than the critical mass Mwr above 
which WR stars form due to the complete removal of the hydrogen envelope. This boundary 
is not well known because of uncertainties in the mass-loss rate; for the solar metallicity the 



currently preferred value is in the range of 30 — 40 M & (IHeger et al.l 120031 ) . Thus, in the 
single star scenario the progenitor mass should be < 40 M & . 

The merger scenario includes two major stages: (i) conservative evolution of a close 
massive binary that ends up with the formation of the neutron star (black hole) in pair with 
the massive normal star, and (ii) spiral-in of the neutron star (black hole) in the massive 
star companion with subsequent me rging with the helium core. The neutron star in the 



process of merging grows black hole (jChevalierlll993l ) so that the final stage of pre-SN is the 



central black hole surrounded with the dense helium accretion disk embedded in the massive 
hydrogen envelope. 

As an example, the original close binary might consist of M\ = 25 M & and M 2 = 20 M & 
stars. After the conservative mass transfer the primary ends up as 5 — 8 M & He star in 
pair with the secondary 37 — 40 M Q star. The subsequent evolution of the helium star leads 
to SN Ib/c or SN lib explosion with the neutron star (black hole) left behind in pair with 
the massive star. This is the stage of massive X-ray binary. After the secondary forms a 
He core, the stellar radius grows, the Roche lobe overflow begins and the common envelope 
forms. The common envelope stage ends up as a black hole surrounded with the helium 
accretion disk within the massive hydrogen envelope. The rapid disk accretion onto black 
hole eventually gives rise to the SN 2009kf event. 

The energetic SNe IIP produced by the hypernova mechanism are extremely rare events. 
This is indicated by the large distance of the unique SN 2009kf. To a first approximation 
the rate ratio of SN 2009kf-like to SN 1998bw-like events can be estimated interpreting their 
distances (Di = 740 Mpc and D 2 = 40 Mpc, respectively) as the closest neighbor distances. 
The closest neighbor distance R is defined by the obvious equation (l/3)R 3 Qr]Gt = 1, where 
Q is the solid angle of the survey, r\ is the detection efficiency, G is the SN rate per unit 
volume, and t is the total observation time. Assuming that Q, r], and t are comparable for 
the detection of both hypernovae and SN 2009kf-like events, we obtain the rough estimate 
of the rate ratio G\jG 2 ~ (D 2 /Di) 3 ~ 10~ 4 . Note that Malmquist bias cannot change 
this ratio significantly because the absolute magnitudes of these phenomena are comparable. 
We conclude therefore that roughly 0.01% of all hypernovae can retain a massive hydrogen 
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envelope and produce SN 2009kf-like events. 
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Fig. 1. — The calculated V light curve of th e optim al model (solid line) is compared with the 
observations of SN 2009kf (IBotticella et al.l (120101 )) (open circles). The effect of the adopted 
56 Ni mass is illustrated for the cases of the 56 Ni mass equal to 0.0765 M (dashed line) and 
the zero 56 Ni mass (dotted line). 
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Fig. 2. — The calculated bolometric light curve of the optimal model (solid line) overplotted 
on the bolometric data of the black-body fit for SN 2009kf (open circles). 
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Fig. 3. — The calculated photospheric velocity (solid line) in comparison with the velocity 
at the photosphere recovered from the Ha and He I 5876 A line profiles on day 61 (filled 
circle). 
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Fig. 4. — Explosion energy a) and 56 Ni mass b) versus hydro dynamic progenitor mass for 
six type IIP supernovae. 
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Table 1. 


Hydro 


dynamic 


models of typ 


e IIP supernovae 


and their 


prog 


;enitor masses 




SN 


Ro 




E 


M m 


..max 
V Ni 


mill 


M NS M pre -SN M lost 


AfzAMS 


ref 




(Re) 


(M @ ) 


(10 51 erj 


? ) (1O" 2 M ) 


(km s~ 1 ) 


(km s — 1 ) 


(Mq) 


(M Q ) 


(Af©) 


(Af©) 




1987A 


35 


18 


1.5 


7.65 


3000 


600 


1.6 


19.6 


1.7 


19.8-22.8 


1 


1999em 


500 


19 


1.3 


3.60 


660 


700 


1.6 


20.6 


1.6 


21.0-23.4 


2 


2003Z 


229 


14 


0.245 


0.63 


535 


360 


1.4 


15.4 


0.2-0.8 


14.4-17.4 


3 


2004et 


1500 


22.9 


2.3 


6.8 


1000 


300 


1.6 


24.5 


1.4-3.4 


25.0-29.0 


1 


2005cs 


600 


15.9 


0.41 


0.82 


610 


300 


1.1 


17.3 


1.0 


17.6-20.4 


5 


2009kf 


2000 


28.1 


21.5 


40.0 


7700 


410 


4.5 a 


32.6 


3.4 


> 36 b 


6 



Note. — Table \T\ gives mass estimates of type IIP progenitors according to the following relation: AfzAMS = Afjvs or Mbh + M env 
+ M lost . 

a It is the mass of black hole 

b In a single star scenario 
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